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Abstract 

It is shown that as far as the linear diffusion equation meets both time- 
and space- translational invariance, the time dependence of a moment of 
degree a is a polynomial of degree at most equal to a, while all connected 
moments are at most linear functions of time. As a special case, the 
variance is an at most linear function of time. 
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1 Introduction 



Anomalous diffusion has found a special place in both observational and theoret- 
ical sides of an extensive list of disciplines including physics, chemistry, biology, 
economy, engineering, geology, meteorology and astrophysics [1—5]. This phe- 
nomenon is characterized by the long-time behavior of the variance of some 
density behaving like 

Var(X) = (X 2 ) - (X) 2 -^Ct a , t -> oo (1) 

with a 1. For ordinary diffusion a = 1. 

There has been extensive efforts to model the anomalous diffusion, among 
them the continuous time random walks (CTRW) [1-3,6,7]. Interestingly, the 
Fourier-Laplace transform of the CTRW's equation may generate, among other 
forms, the fractional diffusion equation, which has the form of an ordinary 
diffusion equation but with non-integer order of derivatives of space or time [5] . 
In recent years, great interest has been attracted to the study of the fractional 
dynamical equations to model the anomalous diffusion behavior [1,8-10]. 

There have been efforts in other directions to explore in what situations 
one may get the anomalous diffusion behavior. In [11], it was shown that the 
ordinary Fokker-Planck equation, which is local, can give anomalous diffusive 
behavior in situations in which there is a non-zero drift velocity. The cases 
with zero drift velocity are explored in [12], and it is shown that very particular 
forms of temporal and spatial dependent diffusion coefficients can generate the 
anomalous diffusive behavior. It has also been shown that particular non-linear 
diffusion equations, which are like the ordinary ones in form, can produce the 
anomalous diffusion [13]. 

The aim of the present work is to show that systems satisfying the conditions 
of linearity, both space and time translational invariance, and the property that 
the state of the system at one time determines uniquely the state of the system 
at later times, do not show anomalous diffusion. Of course it does not mean 
that the lack of any of the above conditions necessarily results in anomalous 
diffusion. Perhaps the most important implication of this result is that it is 
not the appearance of fractional derivatives, but some inhomogeneity hidden in 
the definitions of some versions of such derivatives which result in anomalous 
diffusion. In fact, there are versions of fractional derivative that commutate 
with translation operators (an explicit example will be given), and it will be 
shown that evolutions based on those do not result in anomalous diffusion. 

The scheme of the paper is the following. In section 2 the relation of time 
translation symmetry and fractional time derivative is studied. Section 3 is 
on the solution of fractional time derivative differential equations. In section 
4 the relation of moments corresponding to a distribution, and the Fourier 
transform of that distribution are discussed. In section 5 the time evolution of 
a distribution is discussed which satisfies a linear equation enjoying time and 
space translation symmetry, and the property that the state of the system at 
one time determines uniquely the state of the system at later times. It is shown 
that the connected moments corresponding to such distributions are at most 
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linear functions of time. In section 6, an example is studied which lacks time 
translation invariance and does exhibit anomalous diffusion. The origin of this 
anomalous diffusion is then traced back (by explicit calculation) to the lack of 
time translation invariance. Section 7 is devoted to the concluding remarks. 

2 Time translation symmetry and the fractional 
time derivative 

There are different versions of fractional derivative. Demanding that the frac- 
tional (time) derivative be commuting with time translation, however, makes 
the fractional derivative almost unique. By almost it is meant that as long as 
a specific branch of complex power function is chosen. Denote by To (a) the 
operator which translates in time by a: 

[T (a)V] := 1>(t-a), (2) 

and by D the time derivative. The generator of the group of time translations 
is (-D ): 

T (o) =cxp(-aD ). (3) 

Any operator which commutes with all time translations (in fact with all time 
translations with the translation parameter in an interval containing more than 
one point) commutes with the time derivative as well. The converse is also 
obviously true. So the lack of time translation invariance of an operator (in- 
cluding a fractional time derivative), is equivalent to the non-commutativity of 
that operator with the time derivative. As the time derivative is not degenerate, 
one then concludes that any such operator is a function of the time derivative. 
Denoting such an operator with O, there is a function / such that 

O = /(Do). (4) 

The eigenvectors of Do are exponentials. Denoting by <?(.) a function the value 
of which for the variable t is g(t), one has 

D [exp(s .)] = s [cxp(s .)], (5) 

so, 

0[exp( S .)]=/(s)[exp(*.)]- (6) 

To define the fractional time derivative, one simply takes / to be a power func- 
tion. But then, a power function with a fractional exponent is not unique. To 
characterize it completely, one should specify, in addition to the value of the 
exponent, the branch cut in the complex plain. So is the case for the fractional 
derivative. 

An example of a fractional time derivative defined in the literature, which is 
not commuting with time translation is the Caputo fractional derivative ( ^Do p) 
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[14]: 

( m) ■= f^ITp) f dT (* - ^ (n) to, n - 1 < /3 < n. (7) 

It is seen that 

T a ?Do,3 = b+aDO/3 T a . (8) 

Another example of a fractional time derivative widely used in the literature, 
which is not commuting with time translation is the Riemann-Liouville frac- 
tional derivative ( R b L Do/3) [15]: 

(^ Do ^ )(t):= ____ ^ dr^-r)"-^- 1 ^), n-l</3<n. 

(9) 

It is again seen that 

T a R b L D0£ = fc+aDo/3 T Q . (10) 

One example of a fractional time derivative defined in the literature, which 
is commuting with time translation is the Weyl fractional derivative ( w Do^): 

< " D « *«" - r[„r + i - rOT ] HSWT i * < T - 



d [/3] + 



1 (-oo 

/ drr™- r («^ m »(T + t), (11) 
Jo 



r[m + 1 - r(/3)] df[«+ 1 7 

where m is any nonnegative integer (the right hand side does not depend on it, 
as long as the integrals are convergent), [/?] is the largest integer not exceeding 
(3, and r(/3) is the fractional part of (i: 

r(/3) :=/?-[/?]. (12) 

Direct calculation shows that if s is a complex number of nonpositive real part, 
then 

w Do p [exp(s .)] = exp[-i n r(/3)] s' 3 [exp(s .)] , (13) 

so 

w Do^cxp[-i7rr(/3)]D^. (14) 

where the branch cut for the power function of exponent (3 has been taken in 
the right half plane, for example on the positive real semi axis. 

3 Solution to fractional time differential equa- 
tions 

Consider the equation 

[/(D o )]V = 0. (15) 
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A solution to this equation is 



^) = £d S ^M, (16) 

where C is an arbitrary contour in the complex plane, and g is an analytic 
function. It is seen that with tp like (16), the left hand side of (15) is an integral 
over C of the product [g(s) exp(si)]. If g is analytic, then the product is an 
analytic function of s and the left hand side of (15) vanishes. 

If / is a polynomial of degree n, then the singularities of the integrand in the 
right hand side of (16) consist of only poles, which are the zeros of /. In that 
case only the values of g and some of its derivatives at those zeros of / which 
are encircled by C enter the integral in the right hand side of (16). To be more 
specific, corresponding to a zero Sj only the derivatives of g at Sj of order less 
than dj enter, where dj is the degeneracy of Sj. The most general solution to 
(15) is then 

k dj 

^^^^"expfet), (17) 
j=i i=i 

where djt's are arbitrary. One of course has 

k 

£>=n. (18) 

i=i 

So the most general solution contains exactly n arbitrary constants. 

If / is not a polynomial, then other sorts of singularities (branch cuts and 
fundamental singularities) may arise as well. It may happen then, that the 
general solution to (15) contains an infinite number of arbitrary constants. 

Now consider equation (15) together with an initial value condition 

V>(*o)=Vo, (19) 

and the boundary condition at infinity 

\m\<M, t>T, (20) 

where T and M are some constants. To achieve the boundary condition, one 
would use a contour which does not enter the complex right half plane. The 
solution would then be unique, if and only if the only singularity of the integrand 
in the right hand side of (16) in the complex left half plane is a simple pole. 
That is the case if in the complex left half plane, / has only poles (no other 
singularities) and just one simple zero. 



4 



4 Moments corresponding to a distribution 



Consider a density ip denned on the space. The moments corresponding to such 
a distribution are expectation values of monomials of space coordinates: 



M(a) := (J 








J dVip{r) 


1 -l 



J dVip{r) Y[(r j ) c 



(21) 



where (r 1 , . . . ) are the coordinates of r. Taking the space to be R™, these can 
be expressed in terms of the Fourier transform of the density. Denoting the 
Fourier transform of tp by (Ftp), 



= J 



(JVO(fc) := / dV exp(-ife-r)V(r) 



one has 



M(a) = [(^)(0)]" 1 



n(^) ™ 



(0). 



(22) 



(23) 



The connected moments are defined as the derivatives of the logarithm of the 
Fourier transform of the density: 



d 



M c ( a ):=\]J(i—) MW)] ( (°) 



(24) 



Examples of the connected moments are the expectation values of coordinates: 



dk 



(25) 



and correlations of coordinates: 



<(r* - (r j )) (r l {r 1 ))) = { (i ^) (i ^) (0), (26) 



from which one constructs the variance: 



Var(r) := <(r - <r» • (r - (r»> , 

= -WN^)]}(o). 



(27) 



5 Time evolution of distributions 

Consider a distribution ijj satisfying the equation 

[/(D ,D)]V = 0, 



(28) 
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where D is differentiation with respect to r (the space variable). Space (time) 
translational invariance means that / does not depend on the space (time) 
variable. It is assumed that this is the case. Note that such operators include 
any space operator subject to only translational invariance, including fractional 
space derivatives with the effect that the space Fourier transform is multiplied 
by a function of the Fourier variable. 
The Fourier transform of (28) reads 



[/(D ,ifc)](^) = 0. 



(29) 



The criterion that the solution to this which does not blow up at t — > oo be 
unique (knowing only one initial condition) is that in the complex left hand 
plane /(.,ifc) has only poles and just one simple zero. Denoting that zero by 
E(k), one arrives at 



(7V)(*,*0 = [(^)(0,fe)] exp[t £(*)]• 
Then, using (24) it turns out that 



d 



[M c (a)](t) = \ J] ( i W{ ) Mm\ )% 0) + 1 



= [Mc(a)](0)+t 



n 



E 



(0). 



(30) 

(0), 
(31) 



This shows that the connected moments are at most linear functions of time. 
Specifically, the variance changes linearly with time. These conclusions are valid 
unless the corresponding derivatives of E at k = either are zero or blow up. 
The moment would be constant in the former case, and would be infinite in the 
latter. So no fractional power time dependence would arise. 

A similar conclusion holds for general (not necessarily connected) moments. 
Using (23) and (30), it is seen that 



[M(a)](t) = [M(a)](0)+P a (t), 
where P a is a polynomial of degree at most deg(ct) in t, with 

deg(a) :=Y, a J- 



(32) 



(33) 



Exceptions are again when the derivatives of E at k = blow up, but even in 
that case one does not encounter a fractional power law evolution either. 

As an example of a linear fractional time derivative equation, which is in- 
variant under time and space translations, consider 



(-l) m w D p ^ + D-D^ = 0, 



(34) 
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where the branch cut for the power function with the exponent /?, is the positive 
real semi axis. Using (13), one has 

exp[-i7rr(/3)] [E(k)f - k ■ k = 0, (35) 

or 

[exp(-in) E(k)f = k ■ k. (36) 

Noting that the angle of E(k) is between and (2 it), it is seen that one solution 
of this equation is 

E{k) = -(fc-fc) 1 /' 3 . (37) 

The solution is unique, provided (3 is less than 2. Assuming that's the case, one 
arrives at 

(JV0(t, k) = [(^)(0, fe)] exp[-t (k ■ kf/% (38) 

which shows that there is no power law time dependence with fractional ex- 
ponents, for connected moments, although such moments could be constant or 
blow up, when (5 is not equal to one. 

6 Lack of time translation symmetry, and non- 
constant diffusion rate 

As an example of a fractional time derivative evolution which results in a non- 
constant diffusion rate, consider 

£d 0/3 V> = D-DV, (39) 

where 

/3 =: 1 - e, (40) 
and e is a small positive number. Using (7), one arrives at 

1 

(D tf)(i) + e 7 (DoVO(t) 



(£ Do ^ m,r) 



f(D m,r)+ [ dr (t-r) £ (D^)(r,r) 
Jo 



(D V)(0) Int+ / dr ln(i - r) (D^)(r) 
Jo 



•o(e), (41) 



where 7 is the Euler constant. To obtain a perturbative (in e) solution of (39), 
one considers the three last terms on the right hand side of (41) as sources 
applied to an ordinary diffusion equation. Then the Fourier transform of (39) 
reads 

(D .FV>)(*, fe) + J(t, k) = —k -kFipit, k), (42) 

where 

J(t, k) = Ji(i, fe) + J 2 (i, fe) + o(e), (43) 
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and 

Ji(t,fe) =e 7 (D ^ H )(t), 



J 2 {t,k)=e 



(D ^ H )(0) lni+ / dr ln(i - r) (D^ H )(r) 
Jo 



(44) 



where V' H is the solution to (42) without the source (J) term. It is seen that 
the relation of J\ to ip is time translation invariant, while that of Ji with is 
not. 

The solution to (42) is 



/>oo 

(^)(t, k) = (.FV H )(i, fc ) + / dt ' G(t - t', k) J(t', k), 

Jo 



where 

(FiPi)(t,k) 



Finally, using (27) one arrives at 

[Var(r)] (t) = [Var(r)] (0) + 2 d f - 2 d e 7 ^ + J dt' In i'^ + o 



(45) 



where 

(JV> H )(i,fc) = (JV)(0,fc) exp(-fe-fet), (46) 
and G is the Green's function, satisfying 

(D G)(t, fc) + = k k G(t, k), 

G(0 _ ,fc)=0. (47) 

One arrives at 

G(i,fc) = -0(i) exp(-fc-fci), (48) 
where 0(i) is the Heaviside step-function. One then has, 

(J r ip 1 )(t,k) = e'yk-k f dt' aip[-k ■ k (t - t')} (Tip H )(t' ,k), (49) 
Jo 

(^)(*,fc) = e 7 fc-fc(J"^ H )(0,fc) / dt' exp[-fe ■ fe (i - *')] hit' 

Jo 

- e (fc • fc) 2 [ dt' cxp[-fc • fe (t - t')] 
Jo 

x / drln(i' -T){I-^ H )(t',k), (50) 
Jo 



poo 

/ dt' G(t-t',k)Ji(t',k). (51) 
Jo 



(e), (52) 



where d is the dimension of the space. It is seen that it is only the integral of 
the logarithm in the right hand side, which results in a fractional power of time 
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appearing in the variance: 



[Var(r)](i) = [Var(r)](0) + 



2dt 1 ~ e 




= [Var(r)](0) + 



(l-e)r(l-e) 
2dtP 



(53) 



The last line is in fact exact, as can be obtained by the exact solution of (39). 
The integral of the logarithm in the right hand side of (52) comes from J2, which 
has a time translation noninvariant relation to tp. So the origin of the fractional 
power of time in the variance is lack of the time translation invariance of the 
evolution equation. 

7 Concluding remarks 

It was shown that a system satisfying the conditions of linearity, both space 
and time translational invariance, and the property that the state of the system 
at one time determines uniquely the state of the system at later times, does 
not show anomalous diffusion. If these criteria are met, then the time depen- 
dence of the connected moments are at most linear functions of time, while the 
time dependence of a general moment of order a is a polynomial of degree not 
larger than a. This result holds regardless of the specific operators entering the 
evolution of system, even if they contain fractional time or space derivatives, 
provided of course the above criteria are met. Two examples of fractional time 
derivative evolutions were also explicitly studied. One of them (the Weyl deriva- 
tive) satisfies the above conditions, hence does not exhibit anomalous diffusion. 
The other (the Caputo) does not satisfy the above conditions, and does exhibit 
anomalous diffusion. In the latter case, the anomalous behavior was traced back 
to the term in the Caputo fractional derivative which violates time translational 
invariance. 

Acknowledgement: This work was partially supported by the Research Coun- 
cil of the Alzahra University. 
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